Time Reversal Symmetry Breaking and Spontaneous Currents in 
s-Wave / Normal Metal / d-Wave Superconductor Sandwiches 



Andrea Huck°, Anne van Otter lo^'*'^, and Manfred Sigrisf^ 
a) I. Physikalisches Institut, Universitdt Hamburg, D-20355 Hamburg, FRG 

b) Physics Department, University of California, Davis, CA 95616, USA 

c) Theoretische Physik, ETH-Honggerberg, CH-8093 Ziirich, Switzerland 

We study the physical properties of an s-wave - normal metal - d-wave junction in terms of 
the Andreev bound states in the normal metal layer. The phase dependence of bound states with 
different orientations leads to superconducting states with broken time reversal symmetry for generic 
orientations of the d-wave superconductor crystal. The occurrence of such a state and the associated 
spontaneous supercurrent along the junction is analyzed in the framework of Ginzburg-Landau 
theory and by the solution of the Bogoliubov-de Gennes equations. 
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During the last few years the order parameter (OP) 
symmetry has been one of the intensively debated issues 
in the field of high-temperature superconductivity. A 
growing number of experiments leaves little doubt that 
the basic symmetry of the Cooper pairs has d^2_j^2-wave 
character in many of the high-temperature superconduc- 
tors (HTSC) 0:^. The unconventional symmetry of the 
OP has important implications for the Josephson effect. 
For d-wave superconductors the Josephson coupling is 
subject to an additional phase dependence caused by 
the internal phase structure of the Cooper pair wave 
function. The phase properties of the Josephson effect 
have been discussed within the framework of the gener- 
alized Ginzburg-Landau (GL) |^ as well as the tunneling 
Hamiltonian approach Q . It was found that the current- 
phase relation depends on the mutual orientation of the 
two coupled superconductors and their interface. This 
property is the basis of all the phase sensitive experiments 
probing the OP symmetry. In particular, it is possible to 
create multiply connected d-wave superconductors which 
generate half-integer flux quanta as observed in experi- 
ments §. 

Various interesting phenomena occur in 45°-interfaces 
of dj;2_y2-vBNe superconductors, where one of the nodes 
of the pair wave function lies parallel to the interface nor- 
mal vector (Fig. 1). For an interface to a normal metal 
or an insulator a bound state appears at zero energy giv- 
ing rise to a zero-bias anomaly in the J-y-characteristics 
of quasiparticle tunneling It was also shown that 

in such an interface to an s-wave superconductor the en- 
ergy minimum corresponds to a Josephson phase different 
from or TT Based on Ginzburg-Landau theories it 
was suggested that this is connected with a local break- 
down of time reversal symmetry T The s-wave 
and d-wave OP can form a complex combination, a so- 
called (s-l-zd)-state, close to this 45°-junction. This leads 
to a phase difference of +tt/2 or —tt/2 across the inter- 
face, which corresponds to two degenerate states JlO| , pTt . 
It can be seen from the GL formulation that under this 
condition a spontaneous current flows parallel to the in- 



terface which produces a local field distribution 
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FIG. 1. Schematic view of the SND-junction. The angle a 
denotes the orientation of the d-wave superconductor (crystal 
axis a and b) and (3 the momentum direction of the bound 
state. The currents generated by the bound states tend to 
cancel in the direction perpendicular to the interface, whereas 
they add parallel to the interface to a spontaneous current. 

In this paper we consider a 45°-interface with a nor- 
mal metal between the d-wave and the s-wave supercon- 
ductor, a device which we call the SND-junction [ p^ . 
Also for this configuration a T-violating state appears 
and generates a supercurrent mainly in the region of 
the normal metal. It is our goal to demonstrate that 
this current has a simple and intuitive interpretation in 
terms of subgap Andreev bound states in the sandwiched 
normal metal layer. Let us first outline the basic idea 
for the situation shown in Fig. 1 where a — -k/A and 
the c-axis is parallel to the interface. In terms of the 
phase difference ^p — Lpd — the Josephson current 
carried by a bound state with a specific orientation (3 
can be expanded in harmonics of the phase difference 
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If as //3(<p) = Ii{l3) sin((p) + hifi) sin(2(p) + • • •. In the 
geometry considered, each bound state with orientation 
< /3 < 7r/4 that sees the "+" lobe with phase ipd, has a 
mirror bound state with orientation —(3 that sees the "-" 
lobe of the d^2_y2-^a.vc wave function with phase ^pd + tt. 
As a result, in the total current perpendicular to the in- 
terface, all odd harmonics cancel, and the Josephson cou- 
pling is reduced. The leading term is ^ sin(2(/3) p^ ] 
and the stable ground state with /j^ = is at = ±7r/2 
and, thus, breaks time reversal symmetry. The Joseph- 
son current parallel to the interface, however, has con- 
tributions from the odd harmonics and to leading order 
/|| ^ sin((^). Remarkably, this parallel contribution is 
nonzero in the ground state and constitutes a sponta- 
neous current. 

Let us first consider this property of the SND junction 
on a phenomenological level by means of GL theory. We 
describe the superconducting state by two OP's, rjs (s- 
wave) and rjd (d-wave), which correspond to the local 
pairing amplitudes. The corresponding GL free energy 
T has the general form. 
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+i^{{u,Vsnu^Vd) - {nyfjsTiUyTjd) + c.c.)] , (1) 

where /o is a free energy density, Tcs and Ted are the 
transition temperatures of rjs and rid, respectively, and 
/3s, d, 7i,2, ^s,d, and ^ are real coefficients {^s,d corresponds 
to the zero-temperature coherence length). These coef- 
ficients and the transition temperatures are in general 
different in the three regions of the SND-junction. We 
use n = V — (27ri/$o)A, with vector potential A and 
flux quantum $o = hc/2e. To study the properties of the 
SND-junction we minimize this free energy with respect 
to rjs^d and A. Assuming homogeneity along the interface 
the problem reduces to one spatial dimension which cor- 
responds to the [l,l,0]-direction in the coordinates used 
in J- (x = a and y = b). We call this direction x' and the 
perpendicular ones y' and z. 

We solve the complete set of GL equations numerically, 
for the case in which the coefficients in are identical 
for both OP's and throughout the system. The tran- 
sition temperatures are only different from zero in the 
corresponding superconducting regions. We assume the 
interfaces between the different layers to be completely 
transparent, i.e. the OP's are continuous and have a 
continuous derivative. For our calculation we choose 
(}s= I3d = 1/2, ^s = = 1 (unit of length), 71 = 4/5, 
72 = 2/5 and ^ = 1. This leads to /o = iJ^/87r, where 
He is the thermodynamic critical field at T = 0. We 
fix $o/2-\/27r i?c^s — 4 which corresponds to the London 
penetration depth A at T = in units of ^g. The result 
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FIG. 2. Spatial dependence of the OP in the SND-junction 
(a = TV /4) based on the GL theory. The parameters of GL free 
energy are given in the text. The temperature is T = Tca,d/'2 
and the width L of the normal metal layer is 2 in units of . 



is shown in Fig. 2 for the OP's and in Fig. 3 for the mag- 
netic field and the supercurrent along the j/'-direction. 

Both OP components penetrate the normal metal layer 
(proximity effect) and coexist there in a combination, 
which for the case a = 7r/4 is entirely determined by the 
mixing terms (72/2)(?7*^?7^ -I- rjlr]*^). Within the weak 
coupling approach which we assume to apply, at least, 
within the normal metal layer, 72 is positive 10. This 
term yields the basic cos(2(^)-dependence of the SND- 
junction free energy. It fixes the phase difference between 
and rjd to ip = ipd — ^fs = ±7r/2 in accordance with 
the argument given above. The mixed state has the T- 
violating s ± id-character in the normal metal. 

The supercurrent density follows from T as J = 
—2cdJ-/dA. We find that the current component Jx' = 
J_L vanishes in the stable junction state and that a spon- 
taneous supercurrent flows parallel to the y'-direction 
and generates a magnetic field distribution Bz in and 
close to the metal layer (Fig. 3). Within the GL- 
formulation the supercurrent Jy' = J|| is caused by the 
spatial variation of the two OP components. 
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where we have omitted the diamagnetic part. Note that 
this part of Jyi has essentially the sin (^-dependence antic- 
ipated above. Under symmetric conditions, Jyi depends 
only weakly on x' inside the normal metal layer as shown 
in Fig. 3. The induced magnetic field is screened perpen- 
dicular to the interface on the scale of the London pene- 
tration depth in the superconducting regions by currents 
flowing in the opposite direction. 

Let us turn now to the microscopic view by considering 
the bound state solutions to the Bogoliubov-de Gennes 
equation in the normal metal layer ndi under the sym- 
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with the Fermi momentum hp = {kx, ky, kz) and k'^+k^- 



\k. 



kp. The integral runs over all transverse momenta 
y\T\kz\ < kp and the sum over all possible numbers 
of multiple Andreev reflections n. The factors take 
the suppression due to thermal decoherence and impurity 
scattering into account, 
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sinh[nii^ /^t] ' 
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where we have introduced the normal metal coherence 
length S.T — hvp /{2'KkBT) in the clean limit, the mean 
free path I, and the effective thickness of the normal 
metal layer Ly^ — Lkp/ {kp — ky — fc^)^/^. 



FIG. 3. Spatial dependence of the spontaneous supercur- 
rent and the magnetic field in the SND junction based on the 
GL theory under the same conditions as in Fig. 2 



metric condition, i.e. the d-wave energy gap in D has the 
form Ad — |A|sign(cos[2(6' — a)]), with the amplitude |A| 
equal to that of the gap of the s-wave superconductor in 
S. We take the Fermi momenta in S, N, and D to be 
equal and the transparency of the interfaces to be high. 
Furthermore, we also neglect the suppression of the en- 
ergy gap near the normal metal and assume the pairing 
interaction to be zero in N. 
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FIG. 4. The ground state phase difi'erence ±(/jo as a func- 
tion of orientation angle a for temperatures corresponding to 
C(T)/L = oo, 2, 1, and 0.5 . 

The total Josephson current is a sum over all possible 
(bound) states near the Fermi energy. If the width of 
the normal metal L is smaller than the thermal coher- 
ence length and the elastic mean free path / in N, the 
Josephson current is given by [f6| , 
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FIG. 5. The junction free energy F as a function of the 
phase for Q=7r/4 and £,t/ L= 10, 2, 1, and for q=0.45 7r/4 (one 
minimum) and 0.757r/4 (two shallow minima) at / L=l. 

The simplest case is that of zero temperature in the 
absence of impurities, so that = Z = cxd and all /„ = 1. 
In this limit the sums over n give sawtooth functions of 
the phase difference, saw[(y9j^^] = [(^^.^ -I- 7r]mod27r. We 
obtain the Josephson currents perpendicular and parallel 
to the junction immediately by angular integration. 



^±>^o[(77 + cos(2a))saw((p) 



-(— — cos(2Q!))saw((^ -I- tt)] 



sin(2Q!)[— saw((^) + saw((/3 + tt)] 



(5) 
(6) 



Here A±^ and denote the perpendicular and parallel 
cross-section of the junction, and Jo = ekp/^mnL). Note 
that the current density Jq is inversely proportional to L, 
as in the GL calculation. The junction free energy F{cp) 
is found by integrating I± with respect to the phase. It 
has two degenerate minima at phase differences (po = 
±[7r/2 — cos(2a)], which correspond to a parallel current 
along the junction /|| = ±5 Jo sin(2a). The ground state 
has (s -I- e*'''"c?)-character in the normal metal layer as 
in the phenomenological treatment, again reflecting T- 
violation. 
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FIG. 6. The parallel (dashed line) and perpendicular (solid 
line) Josephson current densities as a function of phase differ- 
ence for a = 7r/4 and temperatures corresponding to £,t / L = 
10,2,1 (decreasing amplitude). 



For nonzero temperature and in the presence of impu- 
rities, we evaluate /||, and the junction free energy 
F numerically. In Fig. 4 the equilibrium phase difference 
1^0 across the junction is plotted as a function of orienta- 
tion angle a for different temperatures £,t/L in the case 
I = oo. We find that time reversal symmetry is broken 
((ySo 7^ 0, ±7r) only for low enough temperatures, or for 
the orientation angle a exceeding a critical value. For 
a = 7r/4, however, tfQ = ±7r/2 for all T < Tcs,Tcd as 
in the GL treatment. The resulting phase diagram is 
completely consistent with the one found by GL theories 
[0 . The result for the junction free energy F is plotted 
in Fig. 5, and I± and /|| in Fig. 6. Both temperature 
and disorder smear the sharp sawtooth structures found 
at T = in the clean limit in a similar fashion. 

The arbitrary equilibrium phase difference leads to ex- 
perimentally observable effects. T-violating junctions 
can lead to phase windings which are non-integer mul- 
tiples of TT, giving rise to non-standard (not (half-) in- 
teger) flux quantization. Thus, it is possible to create 
devices including T-violating junctions which generate a 
spontaneous arbitrary magnetic flux |^,|lO| . The observa- 
tion of such a deviation from standard flux quantization 
is a clear sign of T-violation. Furthermore, the presence 
of two degenerate equilibrium states allows for hysteresis 
effects {ipo <-^- —ipo)- By applying a current through the 
junction one can switch between the two states. This 
effect corresponds to a phase slip with a fractional flux 
moving along the junction This leads to dissipation 
and the enhancement of microwave absorption as soon 
as the junction enters the T-violating phase. The direct 
observation of the spontaneous currents /|| or the field 
might be difficult, since they average to zero over rather 
small length scales (London penetration depth). Thus a 



probe with high spatial resolution would be needed. 

In summary, we have demonstrated that the An- 
dreev bound states in the normal metal layer of an 
SND-junction are the microscopic realization of lo- 
cal T-violation and provide a clear understanding of 
the spontaneous current found in the phenomenological 
Ginzburg-Landau analysis. This observation allows for 
a more quantitative consideration of this effect, which 
will be important for future experimental investigations. 
The experiments discussed at the end are two among sev- 
eral possibilities to observe this T-violating state of the 
SND-junction. Finally, we like to emphasize that this 
effect is only possible in connection with unconventional 
superconductivity and cannot occur for standard SNS- 
junctions. Therefore, high temperature superconductiv- 
ity provides an exciting new class of Josephson phenom- 
ena. 
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